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Wc compare the quenched hadron spectrum on blocked and unblocked lattices for the Wilson quark action, the 
clover action and the tadpole-improved clover action. The latter gives a spectrum markedly closer to the original 
one, even though the cutoff is ~ 500 Mev. 



A great deal of work has been done both on 
perturbative lattice QCD analysis and numerical 
simulations. Yet, recently it has been shown that 
more attention should be paid to simultaneously 
using both tools to make progress in this field [1]. 
Dramatic gains can be achieved by better lattice 
discretization of continuum operators. The Wil- 
son quark action, with its 0(a) discretization er- 
ror, is the foremost candidate for improvement. 

We study such a systematic improvement here, 
testing it on the quenched hadron spectrum. 
Since we do not know what the 'correct answers' 
are for the hadron masses, we use blocked con- 
figurations and try to reproduce the unblocked 
spectrum. Wc use 8"^ x 16 SU(3) gauge configura- 
tions obtained by 2 blocking steps from a 32^ x 64 
lattice at /3 = 6 [2]. The blocking procedure pre- 
serves the infrared properties of the gauge field. 
On the other hand, the quark action is expected 
to show stronger lattice artifacts, which will be 
the object of two improvements: the 0(a) correc- 
tion [3] known as the 'clover' correction, and the 
'tadpole' improvement, which accounts for large 
tadpole renormaHzation of the Hnk variables [1] . 
This will increase the fill-in of the corresponding 
matrix and therefore the computer time needed 
for matrix-vector multiplication. This is more 
than offset by the possibility of going to a coarser 
lattice, where ma is also larger. In our case, the 
gain from the volume alone is a factor 256. 
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For matrix inversion, we use the method of 
stabilized biconjugate gradients BiCGstab2 [4], 
which we show to be faster than conjugate gradi- 
ent (CG) and conjugate residual (CR(1)). 

1. MATRIX INVERSION 

The standard Hnear solvers in QCD simulations 
have been CG and CR in its CR(1) version [5]. 
CG solves the normal equations so that the condi- 
tion number is squared, resulting in slow conver- 
gence. CR(1) is slow compared to its full version 
and convergence is not guaranteed. 

Alternatively, one can use the biconjugate gra- 
dient (BiCG) algorithm, which constructs two se- 
quences of mutually orthogonal residual vectors. 
This is done by applying the nonsymmetric Lanc- 
zos method to the original matrix and to its ad- 
joint. The method does not imply minimization 
of the residual norm ||r|| = ||6 — Ax\\, therefore 
leading to a very irregular convergence history of 
||r||. The newer method of BiCGstab [7] stabi- 
lizes the residual polynomial of BiCG by apply- 
ing to it another polynomial, i.e. VBiCGstab = 
S{A)rBiCG, where S{A} is constructed in factors 
(1 — xi^)(l — X^A).... At each step k, Xk is 
chosen so as to minimize the \\rk\\ of BiCGstab. 
BiCGstab2 does this minimization in the two- 
dimensional complex space of x<;-i and Xk every 
other step, producing faster convergence. In Fig- 
ure 1 it is shown that BiCGstab2 always beats 
CR(1) and CG by a factor ~ 1.3 in the weak and 
strong coupling limits, respectively [6]. 
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Figure 1. The convergence behaviour of different 

matrix inversion algorithms for tiie Wilson quark 
action on an 8'^ lattice, k is adjusted so that 

^max / ^min — 200. 

Further improvement can be obtained by block 
methods, which solve a given linear system for s 
different right hand sides. In this case the dimen- 
sion of the space spanned by the Lanczos vectors 
will be expanded by a factor s so that the number 
of iterations can be expected to decrease by the 
same factor. We have implemented block BiCG 
and see clear improvement on a small lattice. 

2. RESULTS OF SPECTROSCOPY 

We have calculated quark propagators and 
hadron masses for 55 blocked configurations. 
Mass errors are estimated by jack-knife. Using 
the Wilson action as the quark action, the p 
mass was extrapolated through a quadratic fit to 
0.19(3) in units of the unblocked lattice (a'). (To 
rescale the lattice spacing we have matched the 
string tension {aa^= 0.93) to its value on a 32'' 
lattice at 3 = 6 from Rcf. [9]). This value is far 
from the value 0.34(1) extrapolated from the data 
of [8] on the unblocked lattice using the Wilson 
quark action. The discrepancy comes from dis- 
cretization errors on our coarser grid. Therefore, 
we tried to improve the Wilson action systemati- 
cally. 

The 'clover' 0(a) correction [3] to the above 
action gave us rup = 0.25(1), which shows some 
improvement. 

Recently, it has been shown how important tad- 
pole renormalizations are as one tries to connect 



perturbatively the lattice approach to the contin- 
uum theory [1]. This prescription consists simply 
of dividing each link by uq = {\TrUpiaqY^'^ ■ This 
amounts to a trivial rescaling of the hopping pa- 
rameter K kuq in the Wilson action, leaving 
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unchanged (see Table 1). 



We then tried to combine clover and tadpole 
improvement. In this case the extrapolated p 
mass is 0.30(2), markedly closer to the value 
0.34(1) taken as reference. In Fig. 2 p masses are 
compared for the different quark actions, showing 
the important effect of tadpole renormahzation. 
Complete spectroscopy results are presented in 
Fig. 3 for the doubly improved action. 



1.0 



Wilson 
clover 

clover + tadpole 




Figure 2. The p mass for different quark actions 
on 8^ X 16 blocked configurations. Diamonds are 
the results from Ref. [8]. 

The similarity to the unblocked results of [8] is 
remarkable, given that we are working at a cutoff 
a^^ ~ 500 Mcv (equivalent to /? « 5.1). Note 
that we can reach smaller quark masses than in 
[8] at a small fraction of the cost, including the 
regime mjr/mp < 1/2. 

Another indirect measure of improvement is 
the approach of the critical hopping parameter 
Kc to its free value 1/8. The values of Kc listed 
in Table 1 again show the benefits of the doubly 
improved action. 

Improvement of the action is visible as well in 
the eigenvalue spectra. In Fig. 4 we show the 
eigenvalues of matrices corresponding to differ- 
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Figure 3. Hadron masses for clover-and-tadpole- 
improved Wilson quark action on 8^ x 16 blocked 
configurations. Diamonds, squares and octagons 
stand for p, A'', and A masses of Ref. [8]. 



Quark action 


p mass 




Wilson 


0.19(3) 


0.205 


tadpole 


0.19(3) 


0.163 


clover 


0.25(1) 


0.175 


clover+tadpole 


0.30(2) 


0.117 



Table 1 

Changes in Kc and rup for different quark actions. 

ent formulations of quark actions. They are com- 
pared to the spectrum of free quarks on a 32^ 
lattice. 

The similarity of the low-lying spectrum of the 
doubly improved action to the free one is remark- 
able. This result suggests that the free Wilson 
matrix be used as a preconditioner for the lin- 
ear solver. After fixing to Landau gauge, pre- 
liminary results show that the matrix inversion 
can be Fourier Accelerated at least by a factor of 
two, which will be much greater at higher /3. The 
acceleration was also greater as k was increased, 
suggesting that small quark masses can be simu- 
lated efficiently. 

2.1. Acknowledgements 

We thank the QCD-TAR,0 collaboration for 
generating our blocked configurations, and R. 
Gupta for communicating his data to us. 



clover+tadpole 



free ^Vilson 



0.0 0.4 0.8 1.2 1.6 2.0 



Figure 4. Comparison of spectra corresponding 
to different quark actions for k = 0.1. The back- 
ground field in the first three cases comes from a 
4^ lattice at /3 = 6, and the last spectrum is that 
of free quarks on a 32^ lattice with periodic b.c. 
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